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ABSTRACT 

We present the first study of the X-ray population of the intermediate-age rich open cluster NGC 
6819 using the XMM-Newton Observatory. In the past decade, Chandra X-ray observations have 
shown a relationship between the X-ray population of globular clusters and their internal dynamics 
and encounter frequency. We investigate the role dynamics possibly play in the formation of X-ray 
sources in NGC 6819, and compare our results with known properties of field and globular cluster X- 
ray populations. We implement a multi wavelength approach to studying the X-ray sources, utilizing 
X-ray and UV data from XMM observations along with the wealth of photometry and radial- velocity 
data from the WIYN Open Cluster Study (WOCS) and the CFHT Open Cluster Survey. Within the 
cluster half-light radius we detect 12 X-ray sources down to a luminosity of 10^^ erg s~^ for cluster 
members. The sources include a candidate quiescent low- mass X-ray binary (qLMXB), a candidate 
cataclysmic variable, and two active binary systems. The presence of a qLMXB in an open cluster is 
previously unexpected given the known relationships between luminous X-ray sources and encounter 
frequency in globular clusters, and most likely has a dynamical origin. 

Subject headings: binaries: close — open clusters and associations: individual (NGC 6819) — X-rays: 
binaries 



1. INTRODUCTION 



It is well established that dynamical interactions play 
a large role in the formation of exotic binaries in glob- 
ular clusters. From the early days of X-ray astronomy, 
the overabundance of low-mass X-ray binaries (LMXBs) 
in outburst in globular clusters (Clark 1975 Katz||1975 ) 



pointed to formation scenarios un ique to the dense envi- 



ronments within globular clusters. Verbunt & Hut ( 1987) 
showed that the number of luminous LMXBs known 
at the time was in agreement with formation theories 
based on close stellar encounters. After the launch of 
the Chandra X-ray Observatory^ a large number of qui- 
escent LMXBs were found in globular clusters, and sev- 
eral groups independently established that these quies- 
cent LMXBs, like those in outburst , were fo rmed through 
dynamical interactions ( Gendre et al. 2003;|Heinke et al. 
20Q3[|Pooley et al.||2QQ3| r 

Cataclysmic variables (CVs) have also been predicted 
to have an en hancement in their formation rates in glob- 
ular clusters (IVerbunt fc Hut|ll987 Di Stefano & Rap- 
paportp994 Ivanova et al.||2U06[ ), although this was dil- 
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ficult to verify until Chandra b egan finding CVs in l arge 
numbers in globular clusters. [Pooley fc Hut (2006) es- 
tablished that dynamics do play a role in the observed 
population of globular cluster CVs and that some com- 
ponent of the CV population is primordial in origin. The 
relative numbers of the dynamically formed and primor- 
dial populations is subject to a number of poorly un- 
derstood but important dynamical processes. For exam- 
ple, dense clusters would be expected to form many CVs 
through dynamical encounters. However, similar encoun- 
ters would be expected to destroy the wide progenitor s 
of CVs before they entered the CV phase (iDavies '19971) 
Studies of low- density globular clust ers | Kong etdl. 
2006[ |Bassa et al."2008( |Lu et al.|[2QQ9| |Lan et al. 20] 



have revealed evidence that the X-ray source popula- 
tion (mainly CVs and chromospherically active main- 
sequence binaries) in these clusters is likely largely pri- 
mordial in origin with few, if any, dynamically formed 
sources. An extrapolation of this trend further down in 
mass and density to the realm of rich open clusters would 
predict almost zero X-ray sources, although significant 
populations ar e found in M67 ( ,van den Berg et al.„2004) 
and NGC 188 ( |Gondoin|[2QQ5l ) . 

Dynamics may be playing a larger role in open cluster 
environments than globular cluster studies would pre- 
dict as well. The high binary fraction of blue strag glers 
in NGC 188 and th eir unusual orbit parameters (M ath-| 
ieu & Geller 2009) could be indicative that dynamical 
encounters playea a significant role in their production 



dLeigh fc Sillsl|2011[) . Complete AT-body modeling of the 
open cluster 1V167 also points to dynamics playing a sig- 
nificant role in the formation of blue straggler popula- 
tions ( [Hurley et al.|2005| . Further, one of the binaries in 
NGC D^19 may be the result of an e xchange encounter 
dGosnell et ar||2QQ7i [Hole et al.|[2QQ9| . 

'I'his evidence seems to indicate that dynamical inter- 
actions are playing a role in the production of exotic 
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objects in rich open clusters, but the nature of that role 
is only beginning to be uncovered. In order to gain a 
broader perspective on this, we have undertaken a sys- 
tematic survey with XMM-Newton of eight rich open 
clusters previously unobserved in X-rays to be combined 
with archival analyses of four other clusters. Our im- 
mediate aim is to identify and classify the X-ray sources 
through their X-ray, UV, and optical properties, and our 
ultimate goal is to investigate the role of dynamical in- 
teractions in rich open clusters in the production of lu- 
minous X-ray sources, which are usually close binaries. 

We present here our results on the first cluster of our 
study, NGC 6819. Physical parameters for t his cluster 



have been estimated by nu merous authors (see Kalirai et 
al.|2QQlb ; Hole et al.|2QQ9[ and references therein). Here, 



we adopt an extmctio n of E{B — V) = 0.15 given by 
Bragaglia et al. (2001) and convert t his to an equivalent 



column of Nh = 8.3 x 10^^ cm-^ (iPredehl & Schmitt 



|T995). Recent Kepler data were used to determine a 
well-constrained distanc e of 2.34 zb 0.06 k pc and an age 
between 2 and 2.4 Gyr ( |Basu et al]|2011[ ). We take the 
optical c enter of the cluste r to be 19Ml""18^ +40°ll'4r 
(J2000) dHole et al.|[2QQ9| . 

Our X-ray, U V , and optical observations are presented 
in ^ In ^ we bring together the mult i- wavelength data 
on each X-ray source within the half-light radius of NGC 
6819. Our X-ray flux estimates are given in Q and we 
discuss source classification in ^ Our discussion is in ^ 
and our conclusions are laid out in ^ 

2. OBSERVATIONS AND DATA REDUCTION 

We aim to investigate the X-ray sources of NGC 6819 
as thoroughly as possible, and we therefore implement a 
multi wavelength study of X-ray source characteristics. 
In addition to the XMM X-ray data outlined in §2.1[ 
we also include ultraviolet (UV) and optical photometry 
and spectroscopy. We obtained UV photometry using 
the Optical Monitor aboard X MM , the data and reduc- 
tion of which are outlined in 82.21 We also make use of 



the extensive optical photometry provided by the CFHT 
study of NGC 6819 ( [Kalirai et al.||2001bD and the pho- 
tometry and radial- velocity information from the WIYN 
Open Cluster Surve^f] (WOCS, Mathieu 2000^ , as out- 
lined in 82.31 



2.1. X-Ray Data 

NGC 6819 was observed by the XMM-Newton orbiting 
observatory for 25.0 ks (15:35:07-22:32:22 UT) on 2008 
May 18 as part of our X-ray survey of rich open clusters. 
The European Photon Imaging Camera (EPIC) pn cam- 
era and the two MOS cameras were used in full frame 
mode with the medium aluminum filter to block optical 
light. The telescope boresight location was 19Ml"^18^ 
+40°11'12'' (J2000). The field covered by XMM is shown 
in Figure [T] 

We use XMM Science Analysis Software (SAS) version 
10.0 for our data analysis, following the routines outlined 
in the threads on the SAS help pagj^ Data from all 
three cameras are screened for high-energy background 
flaring, although we do not find flaring during the obser- 

9 This is WIYN Open Cluster Study paper XLVII. 

see http://xmm.esa.int/sas/current/documentation/threads/ 




Fig. 1. RGB Digitized Sky Survey image of NGC 6819 overlaid 
with the XMM field of view. The different camera fields of view 
are indicated in the figure. The 3.^3 half-light radius of NGC 6819 
is shown with the cyan dashed circle. The image is 39^ on a side. 

vation. Filtered event lists are created using the evse- 
LECT task to include only those events with energies be- 
tween 0.2 keV and 10 keV. We restrict the filtered events 
using PATTERN <= 12 and <= 4 for the MOS and pn 
cameras respectively, only allowing pixel patterns related 
to valid X-ray events that are well calibrated. Images for 
each camera are created from the filtered event lists using 
EVSELECT. Individual tasks within the edetect_chain 
task are run manually to detect sources. The emlde- 
TECT task uses filtered event lists from all three cameras 
simultaneously, resulting in a combined source list with 
position and photon counts. This combined list contains 
69 X-ray sources in our field. We refer the reader to 
the SAS threads for a more detailed description of these 
routines. 

2.2. Ultraviolet Photometry 

UV data were obtained using the XMM Optical Mon- 
itor (OM), a 30 cm telescope mounted next to the X- 
ray mirrors to allow simultaneous X-ray and UV or op- 
tical observations of the central 17' x 17' X-ray field of 
view. NGC 6819 was observed with the UVWl filter 
(220-400 nm) for 5 ks and with the UVM2 filter (200- 
280 nm) for two exposures of 5 ks each. All data were 
taken in imaging mode. The UVM2 images are aligned 
and stacked using the IRAF routine imcombine in or- 
der to reach a fainter limiting magnitude. We implement 
the SAS pipeline routine OMICHAIN for the OM data re- 
duction. Within OMICHAIN, standard aperture photom- 
etry is applied to the detected sources, taking into ac- 
count the point spread function of the instrument. The 
output of OMICHAIN includes a combined source list for 
3223 sources with calibrated data, including instrumental 
magnitudes, positions, and errors. We use the interactive 
photometry routine OMSOURCE to determine UV magni- 
tudes or limits for all X-ray sources, supplementing the 
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pipeline results for coincident UV objects that are not 
initially detected. 

2.3. Optical Photometry and Spectroscopy 

To identify optical counterparts, we utilize optical data 
from t he the CFHT Open Cluster Survey (Kalirai et ah] 
200 laD and th e WIYN Open Cluster Survey (WOCS, 
Mathieu 2000). Each of these surveys have obtained 
deep, mult i band, wide-field imaging of the stellar pop- 
ulations of NGC 6819. 

For the CFHT Survey, [Ka hrai et al.| ( [2001b[ ) used the 
CFH12K instrument to obtain multiple deep and shal- 
low observations in two filters {B and V) over a 42' x 28' 
field of view that extends well beyond the cluster radius. 
The multiple images are stacked together and PSF pho- 
tometry is performed to measure the position and bright- 
ness of each source. The final color-magnitude diagram 
(CMD) reveals a very tight main- sequence with several 
thousand member stars and extends down to V = 24 
(corresponding to a few tenths of a Solar mass). The ob- 
servations also reveal a rich white dwarf cooling sequence 
in the faint blue part of the CMD w hich stretches over 
four magnitudes (Kalirai et al. 2001b, Figure 7). Further 
discussion of the data analysis of t his cluster, and the de- 
rive d parameters, a re described in Kalirai et al. ( |2001b ) 
and lKahrai fc Tosi| ( |2004[ ) . 

WO (J IS has over 10 years of high-precision radial- 
velocity (RV) measurements over a 23' square field of 
view centered on NGC 6819. The WOCS primary tar- 
get sample includes main-sequence, subgiant, giant, and 
blue straggler stars, spanning a magnitude range of 
11 < V < 16.5 with an approximate mass range of 1.1- 
1.6 M0. Using the Hydra MOS fiber- fe d spectrograph on 
the WIYrs[^3.5m telescope. Hole et alj (|2009|) collected 
RV measurements of ~2700 objects, each with a precision 
of 0.4 km s~^, with observations still ongoing. The pri- 
mary WOCS RV sample with V < 15.0 mag is 89% com- 
plete and the secondary sample down to F = 16.5 mag is 
~60% complete, resulting in 480 RV-determined cluster 
members. These data also include orbital solutions for 
velocity variable stars within the completed sample with 
periods of up to ~1000 days, providing valuable infor- 
mation for possible binary optical counterparts. A more 
thorough d iscussion of the WOCS data are described in 
Hole et aT] ([2009 ) , although we have continued observ- 
ing the cluster since then and include additional RV data 
here. A det ailed descript i on of the analysis techniques is 
provided by |Geller et ai7| ( |2008D . 

3. IDENTIFYING MULTI- WAVELENGTH COUNTERPARTS 
AND CLUSTER MEMBERS 

Twelve of the 69 X-ray sources in the field of NGC 
6819 lie within the cluster half-light radius of 3.3 arcmin. 
T he half-light radius is calculated from the star counts 
in [Kalirai et al.'C2001b ). We only expect five or six back- 
ground sources within this region using t he log A^— log 
S relationship from [Giacconi et al.l ( 2001 ) . Position and 
modeling results (discussed in J4| tor sources within the 
diameter of NGC 6819 are listed in Table [T] and infor- 
mation for sources outside the diameter are listed in Ta- 
ble [2] Sources within the cluster half-light radius are 

The WIYN Observatory is a joint facility of the University of 
Wisconsin - Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatories. 




Fig. 2. — Residuals of the corrected XMM OM astrometry after 
matching XMM UV sources to optical counterparts in the WOCS 
catalog. 

numbered XI through X12 according to their total net 
counts in the EPIC camera. Source position errors are 
calculated using the method of the XMM Serendipitous 
Survey by adding the statistical la errors from emlde- 
TECT to the accept ed systematic XMM position error of 
V!0 in quadrature ( [Watson et aL]|2009D . 



3.1. X-ray Source Cross- Correlation 
3.1.1. Astrometry 

Corrections for systematic offsets between XMM and 
optical positions are necessary before completing success- 
ful source cross-correlation. The WOCS NGC 6819 op- 
tical catalog h as an absolute a strometric error of only 
50 to 100 mas (|Hole et al. 2009), so we use its positions 



as a baseline. We match the XMM frame to the opti- 
cal frame by correlating the OM and optical positions 
through an iterative sigma clipping technique. First, we 
overlay WOCS object positions on the UVWl OM image 
and identify 20 bright counterparts by eye. From these 
sources we compute average offsets in right ascension and 
declination. This offset is applied to all 3223 UV source 
positions and is cross-correlated to the optical catalog. 
Sources are considered a match if they fall within the 
OMICHAIN Icr UV position error. Using the new coun- 
terpart list we calculate the average offset and apply an 
updated correction to the original UV source positions. 
This process is continued until no new counterparts are 
found, resulting in a total of 1150 optical-UV counter- 
parts. The final average offsets are — 3V75 ± 0V14 and 
0V65 ± OVll in right ascension and declination, respec- 
tively. This average offset is applied to all UV and X-ray 
source positions. The residuals of our astrometry correc- 
tion have a sigma of 0^17 and are shown in Figure |2j 

3.1.2. Cross- Correlation 

We aim to classify the X-ray sources using their broad- 
band spectral characteristics, which is heavily aided by 
identifying UV and optical counterparts. After applying 
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TABLE 1 

XMM X-RAY SOURCES WITHIN THE HALF-LIGHT RADIUS OF NGC 6819 

















Error ^ 










Reduced 






a 

oourcG 




Position (J2000) 


(circsGc) 


c 


Lx (ergs 


s-i)d 


IVlOQGi 


Qf of icf 1/^ 


C^-Vciliie 


LyorQnients 


XI 


19 


41 


14.08, 


+40 


14 


05.60 


1.07 


816 


(1.05 + 0.01) X 10^2 


kTNSA = 62.8l°:5 


0.88 


1.0 


qLMXB 


X2 


19 


41 


29.19, 


+40 


11 


23.83 


1.07 


727 


(6.8 + 0.5) 


X 10^1 


kTAPEC = 6.2t2J 


0.96 


0.98 


CV 


X3 


19 


41 


25.41, 


+40 


11 


44.63 


1.16 


361 


(3.7 + 0.5) 


X 10^1 


a— i.O_g g 


0.70 


1.0 




X4 


19 


41 


13.30, 


+40 


12 


48.14 


1.21 


292 


(2.5 + 0.3) 


X 10^1 


kTAPEC = 


0.76 


1.0 




X5 


19 


41 


22.47, 


+40 


13 


06.34 


1.27 


267 


(2.1 + 0.4) 


X 10^1 


kTBB = 0.9t°J 


0.82 


1.0 




X6 


19 


41 


22.02, 


+40 


12 


01.18 


1.40 


155 


(1.6 + 0.1) 


X 10^1 


a=3.0+°-^ 


0.66 


1.0 


RS CVn 


XT 


19 


41 


13.86, 


+40 


09 


32.70 


1.55 


114 


(8.4 + 1.9) 


X 10^0 




1.03 


1.0 




X8 


19 


41 


06.34, 


+40 


11 


42.42 


1.79 


113 


(1.1 + 0.1) 


X 10^1 




0.56 


1.0 




X9 


19 


41 


25.89, 


+40 


12 


23.63 


1.84 


111 


(9.3 + 1.3) 


X 10^0 




0.60 


1.0 


Active Binary 


XIO 


19 


41 


11.81, 


+40 


10 


43.01 


1.96 


91 


(9.1 + 4.9) 


X 10^0 




0.54 


1.0 




Xll 


19 


41 


18.48, 


+40 


14 


15.95 


1.89 


88 


(8.7 + 1.2) 


X 10^0 


a=2.8l°:5 


0.60 


1.0 




X12 


19 


41 


26.64, 


+40 


12 


57.70 


1.76 


86 


(1.4 + 0.5) 


X 10^1 


«=1.4l°o:^3 


0.59 


1.0 





^ Numbered according to their total counts in 0.2 10.0 keV band 
^ Total Icr position error 

^ Net EPIC camera counts in 0.2 10.0 keV band 
^ Unabsorbed luminosity between 0.2—10.0 keV 

^ Model parameter listed for sources X1-X5 based on best fit. Parameters correspond to effective temperature for a neutron star atmosphere (kTNSA) in 
eV, APEC plasma temperature (kTAPEc) in keV, thermal blackbody temperature (kTBs) in keV, or powerlaw index {a). 



the astrometry corrections, the EMLDETECT X-ray source 
hst is correlated to the UV, CFHT, and WOCS optical 
source catalogs within the X-ray position error. Optical 
and UV finding charts for sources within the half-light 
radius of the cluster are shown in Figure [3] 

3.2. Identifying Members 

The UV and optical properties of X-ray source coun- 
terparts are listed in Table [Sj Sources without optical 
or UV counterparts are omitted from this table. An op- 
tical or UV object is considered a probable counterpart 
if it is separated by less than la from the X-ray source. 
Optical or UV sources with a separation between la and 
1.5a are considered possible counterparts; their informa- 
tion is shown in the table with italics. We limit our 
possible counterparts to within 1.5a in order to reduce 
the number of possible spurious counterpart identifica- 
tions that may come from the relatively large X-ray po- 
sition errors. While some X-ray sources associated with 
the cluster may be located outside the half-light radius, 
their frequency becomes much smaller than that of un- 
related background sources. For example, out to twice 
the half-light radius (6.6 arcmin), we detect 28 X-ray 
sources, and we expect between 24 and 32 background 
sources us ing the log N - log S relationship of |GiaccQm 
et al. ( 2QQ1| ). For this reason, we restrict our analysis to 
only those sources within the half-light radius. 

We also restrict possible X-ray source cluster member- 
ship to these sources. Sources within the half-light ra- 
dius have optical counterparts with high proper motion 
membership probability and also include the brightest 
sources in the field outside of obvious foreground stars. 
Optical counterparts of sources outside the half-light ra- 
dius for which RV and proper motion data are available 
are shown to be non- members, as can be seen in Table [3j 
For these reasons we consider it unlikely that we miss X- 
ray cluster members by restricting our sample to sources 
within the half-light radius and therefore focus on these 
objects in the following sections. A "true color" X-ray 



image of this region is shown in Figure [4j Sources within 
the half-light radius are numbered XI to X12 according 
to their detected counts in 0.2-10.0 keV from all three 
EPIC cameras. 

4. X-RAY SPECTROSCOPY 

In order to obtain fiuxes and characterize the X-ray 
sources, we extract and fit spectra for all 69 sources, al- 
though more thorough model fitting is only carried out 
on sources X1-X12. Spectral extraction is completed on 
the filtered event list for each camera using the evselect 
task. Some sources fall between CCD chips or outside the 
field of view of either the MOS or pn, and therefore do 
not have a spectrum from each camera. After extraction, 
the spectra for each source are fit simultaneously using 
Sherpa, the da ta modeling and fitting package in CIAO 
4.3 (Freeman et al.l|2001{ [Fruscione et~aL] |20Q6t ^_We 
implem ent JN elder- ivlead simplex optimization j^lder 
Mead|[T965| using Cstat statistics, a modification of the 



Ciiih| ([T979|. For sources within the cluster 



statistics of _ 

half-light radius that have greater than 200 net counts 
in the EPIC cameras (XI, X2, X3, X4, and X5) we 
do model determination between five different XSPEC 
models: th ermal bremsstrahlung, APEC thermal plasma 
( Smith et a l. 2001), blackbody, neutron star atmosphere 
fFavlo v eF^al.||1991 Zavlin et al. "1996), and power law. 
All source models also take into account absorption, in- 
cluding an absorption component using the minimum 
neutral hydrogen column density of Njj = 8.3 x 10^^ 
cm~^, fixed using the optical reddening value of NGC 
6819 and is frozen for model fitting. The APEC model 
abundances are fixed to s olar to match the ne ar solar 
metallicity of NGC 6819 (Bragagha et al. ||200lD . Back- 
ground spectra are extracted from source-iree regions as 
near to the source as possible while attempting to keep 
a similar distance to the readout node, and are modeled 
simultaneously with the source using an absorbed power 
law. 

For the higher count sources the model with the best 
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TABLE 2 

XMM X-RAY SOURCES OUTSIDE THE HALF-LIGHT RADIUS OF NGC 6819 





Source 


Position Err. 
(arcsec) ^ 


Counts ^ 


Fx(xlO-i4) ^ 
(ergs ciD? s"-*- ) 


PL Index ^ 


XMMUJ 


194157.6+401514 


1.02 


3167 


21.4 + 0.5 


3.10 + 0.04 


XMMUJ 


194216.0+401046 


1.04 


1823 


24.0 + 1.3 


2.07 + 0.05 


XMMUJ 


194046.7+402153 


1.07 


1227 


13.8 + 0.5 


2.95 + 0.07 


XMMUJ 


194121.1+400044 


1.15 


504 


3.5 + 0.2 


3.2 + 0.1 


XMMUJ 


194206.7+401441 


1.15 


472 


5.1 + 0.8 


2.9 + 0.1 


XMMUJ 


194200.0+401520 


1.14 


458 


4.1 + 0.3 


2.7 + 0.1 


XMMUJ 


194057.3+402031 


1.21 


397 


3.7 + 0.6 


1.44 + 0.09 


XMMUJ 


194049.0+400615 


1.22 


378 


2.3 + 0.4 


1.8 + 0.1 


XMMUJ 


194100.2+400220 


1.20 


336 


5.1 + 0.7 


1.7 + 0.1 


XMMUJ 


194146.6+401014 


1.10 


334 


5.7 + 0.7 


2.1 + 0.1 


XMMUJ 


194050.1+401650 


1.20 


331 


3.6 + 0.4 


2.1 + 0.1 


XMMUJ 


194212.7+401747 


1.13 


330 


9.8 + 1.3 


1.68 + 0.09 


XMMUJ 


194030.8+400641 


1.30 


312 


3.4 + 0.8 


3.2 + 0.2 


XMMUJ 


194224.4+400929 


1.37 


309 


3.5 + 1.4 


3.5 + 0.1 


XMMUJ 


194104.8+400100 


1.15 


307 


5.7+1.2 


1.2 + 0.1 


XMMUJ 


194209.9+401449 


1.25 


295 


4.8 + 0.9 


1.8 + 0.2 


XMMUJ 


194141.5+400658 


1.30 


292 


4.8 + 1.7 


0.7 + 0.1 


XMMUJ 


194153.0+401145 


1.17 


274 


5.7+1.2 


1.6 + 0.1 


XMMUJ 


194053.4+401556 


1.32 


247 


5.4 + 3.5 


2.9 + 0.5 


XMMUJ 


194145.4+400806 


1.48 


245 


1.5 + 0.6 


2.2 + 0.3 


XMMUJ 


194137.5+400523 


1.14 


218 


5.5 + 1.3 


1.2 + 0.1 


XMMUJ 


194147.0+401347 


1.34 


203 


3.0 + 1.0 


1.3 + 0.2 


XMMUJ 


194115.8+401724 


1.40 


178 


0.8 + 0.1 


3.3 + 0.2 


XMMUJ 


194034.9+401152 


1.60 


178 


2.6 + 2.8 


3.2 + 0.5 


XMMUJ 


194153.5+401531 


1.52 


177 


2.6 + 0.6 


1.9 + 0.2 


XMMUJ 


194105.5+401430 


1.42 


163 


1.5 + 0.6 


1.7 + 0.2 


XMMUJ 


194131.3+402033 


1.53 


161 


1.5 + 0.6 


1.6 + 0.3 


XMMUJ 


194057.1+395959 


1.36 


157 


0.8 + 0.1 


2.8 + 0.2 


XMMUJ 


194106.6+400753 


2.01 


133 


1.1 + 0.3 


2.4 + 0.3 


XMMUJ 


194200.2+400516 


1.89 


132 


0.5 + 0.1 


2.8 + 0.3 


XMMUJ 


194139.5+402141 


2.07 


128 


0.6 + 0.4 


1.9 + 0.4 


XMMUJ 


194013.1+401143 


1.88 


128 


1.3 + 0.7 


2.2 + 0.4 


XMMUJ 


194226.2+400934 


1.74 


122 


0.7 + 0.5 


2.7 + 0.3 


XMMUJ 


194209.0+401113 


1.70 


122 


1.6 + 0.7 


1.9 + 0.3 


XMMUJ 


194037.0+401141 


1.59 


112 


1.1 + 0.2 


3.3 + 0.2 


XMMUJ 


194121.9+402336 


2.93 


111 


0.9 + 0.2 


3.5 + 0.3 


XMMUJ 


194041.5+400745 


1.68 


109 


1.5 + 0.7 


1.6 + 0.3 


XMMUJ 


194155.3+400940 


1.49 


107 


1.2 + 0.3 


2.6 + 0.3 


XMMUJ 


194149.9+400804 


1.65 


102 


1.6 + 2.4 


3.1 + 0.5 


XMMUJ 


194156.2+401351 


1.92 


101 


2.5 + 2.0 


1.6 + 0.5 


XMMUJ 


194056.2+401541 


1.71 


99 


2.6 + 2.8 


1.4 + 0.2 


Vl\/Tl\/rTT T 


io/in/in oj^/ini^^ii 


i.DZ 


yo 




i.D zt U.z 


XMMUJ 


194037.7+400421 


1.99 


96 


9.2 + 3.8 


1.8 + 0.9 


XMMUJ 


194055.4+401242 


1.87 


92 


0.7 + 0.1 


2.6 + 0.3 


XMMUJ 


194202.4+400734 


2.04 


89 


0.8 + 0.8 


2.5 + 0.6 


XMMUJ 


194141.1+401522 


1.88 


86 


1.6 + 1.2 


2.1 + 0.5 


XMMUJ 


194113.0+400629 


2.45 


85 


5.2 + 2.4 


1.0 + 0.3 


XMMUJ 


194054.2+401036 


1.86 


81 


1.6 + 1.4 


1.3 + 0.5 


XMMUJ 


194021.4+401546 


2.38 


80 


0.3 + 0.2 


1.8 + 0.4 


XMMUJ 


194218.7+401544 


1.77 


78 


1.6 + 0.7 


1.1 + 0.4 


XMMUJ 


194143.4+401600 


1.74 


73 


0.3 + 0.2 


1.6 + 0.4 


XMMUJ 


194130.0+400348 


2.04 


66 


0.6 + 0.2 


3.7 + 0.7 


XMMUJ 


194052.9+400400 


2.06 


66 


0.6 + 0.2 


3.1 + 0.4 


XMMUJ 


194205.7+402102 


2.35 


60 


3.9 + 1.5 


1.9 + 0.3 


XMMUJ 


194200.3+402202 


2.16 


53 


0.9 + 0.2 


3.8 + 0.4 


XMMUJ 


194200.0+401303 


2.40 


53 


0.7 + 0.4 


3.1 + 0.8 


XMMUJ 


194101.7+400411 


1.87 


52 


0.7 + 0.4 


2.1 + 0.5 



^ la error calculated by emldetect 

^ Net EPIC camera counts in 0.2-10.0 keV band 

^ Model-determined absorbed flux between 0.2-10.0 keV 



Best fit power law index 
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Optical Finding Charts 




X1 



UV Finding Charts 



X^X5 

#4 



X7 



o 



9 



o 




m 



XI 



X11 



o 



c3 



X12 



Fig. 3. — V-band and UVWl finding charts for X-ra y sources within the half-light radius of NGC 6819. The optical images are from 
CFHT observations of the cluster (Kalirai et al.|2001bb. Each box is 20^0 on a side overlaid with the corresponding X-ray position error 
circle. The greyscale for the optical X6 chart is lighter Tor clarity. The optical chart for source XIO is shown on a DSS image of NGC 6819 
due to its location between chips in the CFHT images. Sources X7 and X9 are shown using the UVM2 image due to their position near 
dead pixels in the UVWl image. 
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TABLE 3 

X-RAY Source Optical and UV Counterpart Properties 



X-ray ID WOCS ID^ B - UVWl UVM2 PM prob.^ RV prob. ^ Opt. dis.^ UV dis. 



Within half-light radius of cluster 



XI 








19.75^ 


X2 




2140 


0.50 


19.69 


X4 


4003 


13.06 


1.41 


15.43 


X6 


3002 


12.76^ 


1.12 


15.44 


X9 


52004 


15.65 


0.844 


17.28 


Outside half-light radius of cluster 






194157.6+401514 


1017 


6.23 


0.17 




194216.0+401046 


264023 


18.36 


0.59 




194200.0+401520 


64018 


16.58 


1.52 




194121.1+400044 








20.40 


194206.7+401441 


204020 


18.14 


1.55 




194141.5+400658 




20.08 


0.52 


19.73 


194115.8+401724 


3012 


13.21 


0.19 


13.90 


194105.5+401430 


99008 


16.98 


1.65 




194106.6+400753 




22.44 


1.42 




194041.5+400745 


234016 


18.69 


0.38 


19.81 


194140.1+400029 


169025 


17.52 


1.60 




194113.0+400629 




21.33 


1.49 




194121.9+402336 




17.56 


0.76 




194141.1+401522 


255012 


19.29 


0.79 




194132.6+400009 


22024 


14.48 


0.63 




194052.9+400400 




21.73 


1.58 





19.93 








1.45 


20.11 






1.24 


0.79 


17.99 





2 


1.57 


1.64 


16.98 


99 


93 


0.98 


1.16 


18.76 


96 




1.32 


1.14 









0.41 










0.86 










0.61 




19.31 








1.27 








0.96 




21.01 






0.84 


0.90 


15.52 








1.08 


1.23 








1.25 










1.77 




19.69 






1.44 


1.58 








1.98 










0.36 










1.96 










1.58 












1.75 





0.81 



" [Hole et al.1 ( [2QQ9t 



^ Values frorn [KaUrai et aL 



( [20Qlbl , 

Proper motion determined cluster membership probability ( I. Platais, private communication) 
Radial- velocity determined cluster membership probability ( Hole et al. |2Q09| ) 



^ Distance between optical and UV counterparts and the X-ray source are given in arcseconds 
^ Italicized information corresponds to possible counterparts separated from the X-ray source by more than the la X-ray position 




Fig. 4. X-ray "RGB" image of the inner field of NGC 6819. 
The red, green, and blue colors correspond to low, medium, and 
high X-ray energy ranges as indicated in the figure. The magenta 
circle shows the 3.^3 half-light radius of the cluster. The 12 sources 
within the half-light radius are numbered according to their total 
counts in the XMM EPIC camera. The image is 12^ on a side. 



fit statistic and physically reasonable results is reported. 
Goodness of fit is calculated using Cstat statistics, which 
is a suitable reduced statistic method for non-binned 
data. We also report the Cstat statistic Q-value which 
is a measure of the probability that one would observe 
the given statistic value if the model is correct and the 
the best-fit parameters are true. For lower count sources 
with less than 200 net counts, where model determina- 
tion is more difficult, we fit only an absorbed power law 
with absorption in order to estimate the X-ray flux. For 
sources within the half-light radius of the cluster we cal- 
culate the unabsorbed X-ray luminosity from the best-fit 
model using a distance of 2.3 kpc from recent Kepler data 
jBasu et aT.||201l| . We refer the reader to the XSPEC 
manua l for more inf ormation on the implementation of 
Cstat ( |Arnaud|[T996l ). 

The Sherpa task CONF is used to calculate confidence 
interval bounds, which varies a given model parameter 
value while allowing other free parameters to float 
to new best-fit values. Uncertainties for model pa- 
rameters are given as the CONF la confidence limits. 
Flux uncertainties are obtained with the Sherpa task 
SAMPLE_ENERGY_FLUX, in which the model parameter 
confidence intervals are randomly sampled 1000 times 
and the flux calculated for each realization. The stan- 
dard deviation of the resultant flux distribution is taken 
as the l<j flux uncertainty. Modeling results are listed in 
Tables [T] and |2] for the sources within and outside of the 
cluster half-light radius, respectively. 
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5. X-RAY SOURCE CLASSIFICATION 

Six of the 12 sources in NGC 6819 — X5, X7, X8, XIO, 
XI 1, and X12 — are both faint and have no optical or 
UV counterparts. We attempt to classify the other six 
sources — XI, X2, X3, X4, X6, and X9 — based on their 
X-ray and counterpart properties. We note that source 
X3 lacks an optical or UV counterpart as well, but we 
still attempt to classify its emission because it is bright 
enough to carry out model determination. Figures [5^ 
and b show the X-ray counterparts on optical and UV 
CMDs, respectively. In both figures, solid diamonds cor- 
respond to probable counterparts whose distance from 
the X-ray source falls within the la X-ray position error. 
Open diamonds correspond to possible counterparts that 
have a separation between la and 1.5a from the X-ray 
position. 

Previous studies show that studying X-ray sources on 
a diagram comparing luminosity to X-ray hardne ss ratio 



is an effective diagnostic for source classification (Pooley 



& Hut 2006). Hardness ratio here is defined as the num- 
ber of soft counts (0.2-2.0 keV) divided by the number 
of hard counts (2.0-6.0 keV) and is essentially an X-ray 
color, so this diagram can be thought of as an X-ray 
color-magnitude diagram. This type of diagnostic sep- 
arates objects whose X-ray emission mechanisms have 
characteristically harder spectra, such as cataclysmic 
variables (CVs), or softer spectra, such as quiescent low- 
mass X-ray binaries (qLMXBs). In Figure [6] we show an 
X-ray CMD for the 12 sources within the half-light radius 
of NGC 6819, shown as red squares. We also show se- 
curely identified objects from glo bular cluster studi es of 
various types as comparison (see P ooley fc Hutl 2006', and 
references therein). Based on the extensive Knowledge 
and confident identifications of X-ray sources in globular 
clusters, this diagram can be separated into population 
I, II, and III objects. The division between population I 
and II objects is defined for Chandra sources in Pooleyj 
& Hut (2006) and has been converted to an XMM color 
here. Although there is some overlap between regions, 
population I objects are dominated by qLMXBs, pop- 
ulation II objects are largely CVs but have some con- 
tribution from the most luminous active binaries, and 
population III objects are low-luminosity systems that 
most likely include contributions from both population 
I and II as well as millisecond pulsars (MSPs) and RS 
CVn or other active binary systems. 
Another diagnostic tool for determining the source of 

X- ray emission is comparing X-ra y to optical luininosit y 
ratios. In Figure [7| adapted from'Verbunt et al.'C2008'), 
we show the ratio of X-ray luminosity to absolute V mag- 
nitude for the NGC 6819 sources. The optical absolute V 
magnitudes correspond to the X-ray optical counterparts 
detailed in Table|3]or are shown as an arrow at the CFHT 
data absolute magnitude limit of V = 12.76 at the dis- 
tance of NGC 6819. The X-ray luminosities for sources 

XI- X12 are given as the 0.5-2.5 keV luminosity. The 
dashed lines separate different X-ray populations. From 
Verbunt et al. (2008) a line of constant X-ray to opti- 
cal luminosity ratio given by log I/o.5-2.5kev(erg s~^) = 
36.2 — 0.4Mv separates qLMXBs (above the top line) 
from CVs, and a line given by log Lo.5-2.5kev(erg s~^) = 
34.0 — 0.4Mv separates CVs from magnetically active 
binaries (below the bottom line). The solid line is an 



approximate upper bound to the X-ray to optical lumi- 
nosity ratio for active binaries in the solar neighborhood 
given by log Lo.5-2.5kev(erg s-^) = 32.3 - 0.27Mv. For 
mo re information o n the basis of these separations please 
see Verbunt et al. |2008) and references therein. 

'Ihrough combining the information from these clas- 
sification methods we conclude that XI is a candidate 
qLMXB, X2 is a candidate CV, and X6 and X9 are RS 
CVn systems. We are unable to classify X3 and X4 with 
these methods. 

5.1. Candidate qLMXB 

In Figure [Slwe show a background-subtracted pn cam- 
era (2-10 keV) light curve for XI (500s bins). The appro- 
priately scaled background light curve is shown in grey. 
The light curve errors are calculated from the source 
counts using Gehrels approximation to confidence limits 
for a Poisson distribution (Geh rels 1986J . The average 
number of counts for XI in the pn camera is 7.9 counts 
per bin, or a rate of 0.016 counts s~-^, and is shown with 
the dashed line. Analyzing the resulting power spectrum 
yields no coherent periodicity in the source. We also 
tested the light curve against a constant model, which 
yields a reduced = 1-2 with 50 degrees of freedom {p- 
value = 0.16), indicating that XI is consistent with being 
constant. However, visual inspection of the light curve in 
Figure [8] reveals a possible feature around ~7500 seconds. 
That bin contains 15 counts; the probability of obtain- 
ing 15 counts per bin from a Poisson distribution with a 
mean of 7.9 counts per bin is only 1.6%, indicating that 
the feature is likely due to stochastic variability. Based 
on the ambiguous statistical evidence, we can say only 
that XI is mostly constant except for a possible mild 
event in which the count rate roughly doubled for a pe- 
riod of -1000 s. 

The X-ray spectrum of source XI is remarkably soft; 
almost all of its flux is below 2.0 keV, which strongly 
suggests that XI is not a CV. Regardless, we explore 
models typical of both CVs ( APEC thermal plasma) and 
qLMXBs (neutron star atmosphere). We find that no 
APEC model gives a reasonable fit to the X-ray spec- 
trum. The best fit was both statistically unacceptable 
with a reduced of 3.02 for a pn spectrum with 15 
counts per bin and a temperature of only kT = 0.3 keV, 
which i s well outside the range of typical C V tempera- 
tures (Byc kling et"ar][2Q10( Fertig et al.'Ml). 

In contrast, a neutron star atmosphere fits quite well. 
We fit the X-ray spectrum for XI with a neutron star at- 
mosphere model using the canonical neutron star values 
of a radius of 10 km and a mass of 1.4 M©, setting the 
model normalization using the known distance to NGC 
6819. It is not necessary to include a hard power law 
component in the model to achieve a reasonable fit. This 
is in agreement with the spectral properties of globular 
cluster qLMXBs which also lack a hard power law com- 
ponent, in contrast to field qLMXBs which often require 
an addi tional hard power law for successful spectral fit- 
ting (Heinke et al. 2003). Fitting XI with a neutron 
star atmosphere model with an absorption component 
(as described in Q results in an effective temperature of 
62.8±0.5 eV, which is within the typ ically accepted tem - 
perature values for qLMXB systems ( Heinke et al. poo's ). 
The pn spectrum and best-fit model with residuals are 
shown in Figure [9] The data are binned for plotting pur- 
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Fig. 5. — NGC 6819 optical CMD (a) based on the photometry from'Kahrai et al.« ('2001b') and UV CMD (b) based on sources that are 
detected in both UVWl and UVM2 filters. Solid red diamonds correspond to X-ray sources within the diameter of NGC 6819 that have 
optical counterparts that fall within the X-ray position error. Open red diamonds show possible optical counterparts that have a separation 
between Icr and 1.5cr from the X-ray position. Note that source XI has a possible UV counterpart but does not have an optical counterpart. 
Blue diamonds show X-ray counterparts for sources outside the diameter of NGC 6819, many of which have optical counterparts but no 
UV counterparts. 



poses only and are shown with errors calculated using 
data variance statistics. Residuals, calculated as the 
difference of the data and the model divided by the un- 
certainty, are plotted with error bars of unity to show 
the deviation in sigma between the model and the data. 
We note there is an unusual absorption component at 
^0.5 keV, simil ar to what was seen in Aquila X-1 in qui- 
escence by Rutledge et al. (2002). Examination of the 
MOS spectra for XI also show the same absorption. The 
absorption feature cannot be explained by a neutron star 
atmosphere model alone, but perhaps may be indicative 
of metals present on t he surface of the neutron star due 
to active accretion (Rutledge et al. 2002). However, fur- 
ther observations are necessary to determine whether the 
source of the absorption is inherent to the source and any 
possible physical significance it may have. 

XI is well separated from the other X-ray objects in 
Figure |6j lying within the area for population I, X-ray 
transient and qLMXB, objects. There is a possible UV 
counterpart for XI located 1V45 (1.3<j) from the X-ray 
position. The UV source is located in the region blue- 
ward of the main sequence in Figure [5]3, also known as 
the UV-excess region. This region is prima rily domi- 



na ted by objects with accretion disk emission (Dieball et 
al.||20 07), and is a plausible location in a UV CMD tor a 
qLlVlXB counterpart. Although this is also the location 
in a UV CMD where CV populations are found, we note 
that the lack of coherent periodicity as well as the ab- 
sence of a hard component to the X-ray spectrum makes 
it highly unlikely that XI is an intermediate polar CV, 
so we tentatively classify XI as a candidate qLMXB. 



The nearest optical source has V = 16.4, B — V = l.b 
and is located 1^.^56 away from XI (Figure |3|. The UV 
counterpart and nearby optical source are separated by 
0^34, which is a ^2cr separation given the UV and optical 
source position errors, making the association of the two 
rather unlikely. In addition, the optical source is far from 
the main sequence of NGC 6819, indicating that it is not 
a cluster member, and its red color is difficult to reconcile 
with UV emission. The most plausible scenario is that 
the UV source is associated with XI (separated by 1.3cr) 
and the optical source is an unrelated field star. This 
nearby bright star makes detection of the true optical 
counterpart difficult, and it is likely that the true optical 
counterpart to XI is hidden in the glare from this bright 
source. Although we do not detect an optical counter- 
part, the X-ray to optical luminosity limit for XI also 
shows it to be in the range of X-ray transient objects 
(Figure [t]). We note that the lack of an optical coun- 
terpart IS not surprising as qLMX Bs can be extremely 
optically faint ( Heinke et al.||2003 ), with the added com- 
plication of a nearby bright source making detection of 
an optical counterpart very difficult. 

The X-ray spectrum. X-ray color, and limit of the X- 
ray to optical luminosity ratio all point to XI as a can- 
didate quiescent low mass X-ray binary. In addition, 
the absence of a hard component to the X-ray spectrum 
makes it highly unlikely that XI is a dwarf nova or inter- 
mediate polar CV. With this candidate qLMXB classifi- 
cation, XI may be the first system of its kind found in 
an open cluster environment. Higher spatial resolution 
X-ray and optical data are needed to confirm the coun- 
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Fig. 6. — X-ray CMD for sources within the half-hght radius of 
NGC 6819, shown as fihed red squares. The x-axis shows the hard- 
ness ratio, or X-ray color, of the objects. Soft counts are calculated 
from 0.2-2.0 keV and hard counts are calculated from 2.0-6.0 keV. 
Securely identified globular cluster sources are shown as brown 
stars for LMXBs, purple circles for possible LMXBs, blue triangles 
for CVs, inverted gree n triangles fo r pulsars, and orange diamonds 
for active binaries (see |Pooley Hu t 2006 and references therein). 
From these identifications this diagram is roughly separated into 
three different X-ray populations of qLMXBs (population I), CVs 
(population II), and low-luminosity systems (population III) that 
has contributions from I and II as well as most active binary sys- 
tems. The lines separating the populations are defined in Pooley &;| 
|Hut| (^006 ) for Chandra sources and are converted to XMM values 
here. Errors bars are shown in red for XI, X2, X6, X9, and X12, 
and are representative of other sources in the area. The luminosity 
errors for XI, X2, and X6 are too small to be seen on this plot. 
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Fig. 7. — X-ray to optical luminosity ratio for the X-ray sources 
within the half-light radius of NGC 6819, adapted from Verbunt et 
|al.| ( |20 08). Optical luminosity is given as the absolute V magnitude 
ot tne optical counterpart and X-ray luminosity is given as the 
0.5-2.5 keV luminosity. Sources without optical counterparts are 
shown with arrows at the CFHT absolute magnitude detection 
limit of y = 12.76 at the distance of NGC 6819. The arrows in the 
bottom left, from top to bottom, correspond to sources X8, X5, 
Xll, XIO, X7, and X12. The dashed lines of log(Lx/ivy) = C 
separate various source types. The top line {C = 36.2) separates 
qLMXBs from CVs, and the bottom line {C = 34.0) separates 
CVs from active binaries. The solid line is an approximate upper 
limit to the X-ray to optical luminosity ratio for active binaries in 
the solar neighborhood. Solid diamonds correspond to confident 
optical counterparts, while open diamonds correspond to possible 
counterparts that have a separation of Icr to 1.5cr from the X-ray 
source. 
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Fig. 8. — Background-subtracted light curve for XI from the pn 
camera, with 500 second bins. The errors are calculated using 
Gehrels approximation to confidence limits for a Poisson distribu- 
tion (Gehrels 1986). The appropriately scaled background light 
curve is shown in grey, and the average source rate of 0.016 counts 
is shown with the dashed line. Statistical tests of the light 
curve against a constant model show XI to be mostly non-variable 
except for a possible event at ~7500s. 
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Fig. 9. — Spectrum from the pn camera and best-fit neutron 
star atmosphere model for XI. The data are binned for plotting 
purposes only with 15 counts per bin. The errors are calculated 
using data variance statistics. The residuals shown are the 
difference between the data and model divided by the uncertainty, 
plotted with error bars of unity to show the deviation in sigma 
between the model and the data. We note that there is an unusual 
absorption component at ~0.5 keV that is not explained by the 
neutron star atmosphere model. 

terp art identification and sourc e classification. We note 
that [van den Berg et al. (|2004|) discovered a highly vari- 
able and soft X-ray source in M67 (CX 2). They lack a 
confident classification due to the variable nature of the 
source, but detect system parameters that fall between 
expected values for a bl ack hole qLMXB and a neutron 
star qLMXB. However, |van den Berg et al.| emphasize 
the need for more information before final classification. 

5.2. Candidate CV 

Source X2 has a confident UV source counterpart at 
0V79 (0.74cr) and a possible optical counterpart at a dis- 
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Fig. 10. — Spectrum from the pn camera and best-fit APEC 
model for X2. The data are binned to 15 counts per bin for plot- 
ting purposes only. Errors are calculated using data variance 
statistics. The residuals shown are the difference between the data 
and model divided by the uncertainty, plotted with error bars of 
unity to show the deviation in sigma between the model and the 
data. 

tance of 1^/24 (l.lSa). The UV counterpart to X2 is 
located in the blue UV-excess region of Figure [sJd, an 
area dominat ed by objects with a ccretion disk emission 
such as CVs (iDieball et al.||2007|). The possible optical 
counterpart is located in the "gap" region between the 
main-sequence and white dwarf cooling sequence (Fig- 
ure pa), an area also dominated by CV systems. Based 
on the X-ray color and luminosity of X2 (see Figure [6| 
it lies on the boundary between population I and II od- 
jects, indicating the possibility that this object is a CV 
candidate or qLMXB. Its X-ray to optical luminosity ra- 
tio (Figure [7| is more indicative of a CV than an LMXB, 
therefore we classify X2 as a CV candidate. We fit the X- 
ray spectrum of X2 with an APEC thermal plasma model 
with an absorption component from the neutral hydro- 
gen column density to the cluster, resulting in a best 
fit temperature of kT = 6. 2^^ J keV, in good agreement 
with expected dwarf nov ae temperatures ( Byckling et al.j 
2010 iFertig et al. 2011). The pn spectrum and best-tit 
model with residuals for X2 is shown in Figure p^O} binned 
to 15 counts per bin for plotting purposes only. The error 
and residual calculations are the same as in Figure |9] 

The background-subtracted light curve for X2, with 
500s bins, is shown in Figure 11 The appropriately 
scaled background light curve isshown in grey. Errors 
are calculated using the same method as in Figure 
The average number of counts per bin is 5.8, or a rate 
0.012 counts s~^, shown with the dashed line. We find 
a bin with 13 counts, or a rate of 0.026 counts s~^, as 
can be seen at ~10500s. The probability of finding a bin 
with 13 or more counts from a Poisson distribution with 
a mean of 5.8 is only 0.7%. A chi-square test of the light 
curve for X2 with 1000 second bins against a constant 
model yields a reduced = 1.8 with 50 degrees of free- 
dom (p- value = 0.0005), which is inconsistent with a con- 
stant source. However, analyzing the power sp ectrum for 
X2 does not reveal any coherent periodicity. "Eracleous 
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et al. ( 1991 ) show that CV X-ray spectra are often non- 
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Fig. 11. — Background-subtracted light curve for X2 from the 
pn camera, with 500 second bins. Errors are calculated using the 
same method as in Figure [s] The appropriately scaled background 
light curve is shown in grey, and the average source rate of 0.012 
counts s"-*^ is shown with ttie dashed line. A chi-square test of this 
light curve with 1000 second bins reveals X2 is not consistent with 
a constant source, although analyzing the power spectrum does not 
reveal any coherent periodicities. 

constant but lack a clear periodicity, so we take this light 
curve behavior to be consistent with the classification of 
X2 as a CV candidate. 

5.3. Active Binaries 

X-ray source X6 has optical and UV counterparts at 
distances of 0.98 arcsec (0.7<j) and 1.16 arcsec (0.8<j). 
Source X9 has optical and UV counterparts at distances 
of 1.32 arcsec (0.7cr) and 1.14 arcsec {0.6a). Both sources 
are confident cluster members due to radial- velocity mea- 
surements and/or proper motion membership probabili- 
ties (Hole et al. 2009; I. Platais, private communication). 
Their X-ray to optical luminosity ratios are both indica- 
tive of active binary systems (Figure [t]), and their optical 
and UV color and luminosities are consistent with this 
classification (Figure [5). 

It is already known that the optical counterpart to X6, 
woes ID 3002, may have a dynamical history. WOCS 
3002 is a circular binary with a period of 17.7 days that 
appears in the red-clump of NGC 6819 (as seen in Fig- 
ure [5|i). If the primary is a horizontal branch star, the 
orbital separation of this system is too small to permit 
the primary star to exist at the tip of the red giant phase 
without significant mass transfer. As such, |Hole et al.| 
(2009) suggest that the system is the result of a dynam- 
ical encounter that exchanged the current primary star 
into the binary after its red giant stage, and studies have 
shown that this type of exchange is possible ( Gosnell et| 
al. 2007). Since it has an evolved primary star, we clas- 
sify X6 as an RS CVn system. 

Proper motion studies for the optical counterpart to 
X9 give a membership probability of 96% (I. Platais, pri- 
vate communication). The primary star in this system 
is a rapid rotator, which makes finding an orbital solu- 
tion from the WOCS radial-velocity studies extremely 
difficult. However, its X-ray to optical luminosity ratio 
and X-ray emission properties indicate an active binary 
system that we take to be a cluster member due to its 
proper motion. On the UV CMD (Figure (st), source 
X9 is on the red edge of the region populated by main- 
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sequence st ars and red giants (see Figure 4 in Dieball 
et al. 2007). X9 appears red of the main- sequence and 
is brighter than the equal-mass binary sequence on the 
optical CMD (Figure |5k) . The optical location would 
seem to indicate a possible sub-subgiant binary system. 
Sub-subgiants are known to exist in other open clusters 
as well, and have been shown to be a source of X-ray 
emission (IMa thieu et aH2003j [van den Berg et aH2QQ4l ) . 
We classify X9 as an active binary, noting that it is sini- 
ilar to an RS CVn but further observations are needed 
to verify the nature of the sub-subgiant system. 

5.4. Unknown Objects 

Source X3 is difficult to classify. The X-ray emission 
places X3 close to the intersection of all three population 
types in Figure [6) preventing a confident classification 
based on the X-ray information alone. With no opti- 
cal or UV counterpart classification becomes even more 
difficult. The lack of any counterpart may indicate X3 
is a background source, in which case our luminosity 
measurement is only a lower limit. In addition, both 
blackbody and APEC thermal plasma spectral fits yield 
non-physical temperatures of 0.3 keV and 17 keV, re- 
spectively, when compared to the source hardness ratio. 
Without more information we cannot determine the ori- 
gin of the X-ray emission for X3 other than it is best fit 
with a power law index of LSlg g- 

The closest optical/UV counterpart for X4 is WOCS 
ID 4003, located at a distance of 1^57 (1.3<j) from the 
X-ray position. Radial-velocity studies of WOCS 4003 
show it is not a short-period (P < 3000 days) binary, 
meaning it could be a single star or a long period bi- 
nary. WOCS 4003 has a low probability of being a clus- 
ter member from its radial-velocity, and proper-motion 
studie s indicate 0% probability that it is a member ( [Hole 
et al.| (|2009); I. Platais, private communication). Fhis 
object falls near the intersection of all three populations 
in Figure [6| and in the area dominated by active binaries 
in Figure]7| if WOCS 4003 is the true optical counter- 
part and assuming it is at the distance of NGC 6819. 
Since cluster membership is highly unlikely, the luminos- 
ity given is almost certainly incorrect. With only this 
information, the source of the X-ray emission for X4 is 
still unknown. 



6. DISCUSSION 

Without a previous systematic study of the X-ray 
sources in many open clusters, it is difficult to predict 
the types and number of sources we expect to find. If 
we extrapolate the relation between the number of X- 
ray sources and dynamical frequency in globular clusters 
to the regim e of open clusters, we would expect almost 
no sources (Pooley fc Hut 2006^ see below). Studies of 
two other indi vidual open clusters have shown this not 
to be the case (van den Berg et"aL]|2004| |Gondoin||2005| ) 
but little is known about the formation of such popula- 
tions in open cluster environments. Here we make some 
reasonable assumptions about the formation of some of 
our sources, but definitive conclusions await accurate N- 
body modeling of NGC 6819 and systematic analyses of 
a sample of open cluster X-ray populations. 

The presence of a qLMXB in an open cluster may be 
unexpected, and could be assumed to be dynamical given 



the known relationship between dynamics and qLMXB 
populations in globular clusters (Pooley & Hut 2006). 
Although there are currently studies underway to deter- 
mine t he spatial density of qLMXBs in the field ([Grindlay 



et al. 2005), without such a measurement we are unable 
to determine how many primordial qLMXBs would be 



expe cted in an open cluster using field values. [Pooley 



Hut ( 2006) investigate X-ray sources in globular clusters. 



fitting a two-component model to (number of X-ray 
sources per unit mass) vs. 7 (encounter frequency per 
unit mass) that separates the contribution from primor- 
dial and dynamical populations. The model is applied 
to population I, II, and III objects separately as well as 
to all X-ray objects combined. For population I objects 
(dominated by qLMXBs), the best fit to this relation- 
ship results in 0.41q 55 primordial sources per 10^ Mq. 
Although this is a poorly constrained value, we apply it 
here to roughly determine the probability of findi ng a 
pri i nordial population I object in NGC 6819. Kalir ai et| 



al. (2001b) integrate the mass function of NGC 6819 to 
Una a cluster mass of ~2600 M©. Using this value with 
the upper limit for primordial X-ray sources expected 
per unit mass, we expect ~0.002 primordial population 
I objects in NGC 6819. The low likelihood of finding a 
primordial qLMXB in NGC 6819 given this relationship 
would seem to indicate a dynamical history for XI. 

We apply this comparison with caution, noting that the 
validity of extrapolating globular cluster characteristics 
down to the regime of open clusters is dependent, in part, 
on whether the age difference between globular a nd open 
clusters affe cts the observable X-ray population. |Ivanova| 
et al. (2008") explore the number of qLMXBs in globular 



clusters over time, finding that the number of observable 
qLMXBs from 4 to 11 Gyr remains relatively constant 
for the majority of globular clusters, assuming that all 
qLMXBs that appear before 4 Gyr are primordial in na- 
ture. This interesting result would seem to indicate that 
age does not play a major role in the number of observ- 
able qLMXBs, but the initial assumption that all sources 
prior to 4 Gyr are primordial limits the ability to confi- 
dently extrapolate their population characteristics back 
to the age of NGC 6819. Regardless of the formation 
mechanism, the pre sence of a qLMXB is also unexpected 
given the results of Hurley et al. ( 2005 ) that found with 
their A/'-body model only two neutron stars remaining 
in M67 at 4 Gyr, both of which are single stars. Given 
the comparison tools currently available that all point 
to the unlikelihood of creating a primordial qLMXB in 
open cluster environment, we tentatively interpret XI to 
be a dynamically formed candidate qLMXB. 

Predicting the number of expected CVs in open clus- 
ters is similarly troublesome. Although the CV pop- 
ulations of globular clusters and the field have been 
well studied, it is unknown whether either population 
is strictly applicable to rich open clusters at 2 Gyr. We 
extrapolate the relationship for primordial CVs in globu- 
lar clusters down to the masses of open clusters. |Pooley| 
|fc Hut| ( |2006| find a best fit of 0.27li;9^ primordial pop- 
ulation II objects (dominated by CVs) per 10^ Mq. At 
the upper limit, this results in ~0.004 primordial CVs in 
NGC 6819. A more suitable comparison tool may be the 
field CV density. The ChamPLane survey finds a local 



CV space density of 0.9 



+li X 10- 



pc ^ at 95% confi- 
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dence, whic h is consistent with models of the galactic CV 



population ( |Rogel et a'L]|2QQ8[ ). This value is also consis- 
tent within two standard deviations to the local CV space 



density of|Patterson ( 1998 ) at all scale heights. Given the 
l ocal mas s density of the galactic plane of 0.1 Mq pc~^ 
(Kuijken fc Gilmore^l989), this results in a CV mass den- 



^GVs M-^ 



For the mass of NGC 



sity of 9.0+^^ X 10" 

6819 this results in a prediction of 0.23^0.13 primordial 
CVs in the cluster. The upper limit of this prediction, 
0.62 CVs in NGC 6819, is not i nconsistent wi t h a de tec- 
tion of one CV in the cluster. 



(2005) did 



'Hurle y et al.[ 

create one primordial CV in M67 at a cluster mass of 
only 1400 M©, which may point to primordial CVs be- 
ing more prevalent in open clusters than either of these 
comparisons predict. 

The possible formation scenario for the optical coun- 
terpart of X6, outlined in §5.3[ may point to a dynamical 
history, but requires more observations and modeling to 
definitively determine the formation mechanism. In ad- 
dition, the possible sub-subgiant counterpart to X9 can- 
not be explained by single-star evolutionary theory. Pos- 
sible explanations require under- or over-luminous binary 
components that may be formed through mass trans- 
fer, stellar me rgers, and/or dynamical stellar exchanges 
( Mathieu et al.|2003p . Thus, while not yet definitive, sev- 
eral X-ray sources in NGC 6819 hint at dynamical for- 
mation processes. The presence of dynamically formed 
X-ray sources in open clusters would indicate higher dy- 
namical formation rates than seen in globular clusters. 

Primordial X-ray sources in open clusters must also 
follow different trends than those seen in globular 
clusters. Studies show that the X-ray sources of low 
density globular clusters scale more strongly with cluster 
mass than encounter frequency, which would i ndicate a 



populati on dominated b y primordial sou rces ( [Bassa et 
[al.(|2008| |Lu et a l.|[2QQ9{ |L an et a l." 

lactor for the number of primordis 



2010|. The scaling 

primordial X-ray sources per 
unit mass, however, is low and predicts no close binary 
X-ray sources in the mass regime of rich open clusters. 
It is possible that the low density environment of 
open clusters allows many more close binary products 
formed to persist, rather than be destroyed by the 
dense environment of globular cluster centers (|Davies| 
[1997). We encourage further study into the theoretical 
modeling of creating close X-ray binaries in open cluster 
environments so that these formation pathways can be 
better understood. 



7. SUMMARY 

We conduct the first X-ray study of the rich open clus- 
ter NGC 6819, the first of eight open clusters in a sur- 
vey of X-ray populations in open clusters using XMM- 
Newton. Incorporating the use of Optical Monitor UV 
data and extensive optical photometry and spectroscopy 
from the CFHT Open Cluster Survey and WOCS, we im- 
plement a multi wavelength approach to classifying the 
X-ray population of NGC 6819. We find 12 X-ray sources 
within the half-light radius of NGC 6819 to our luminos- 
ity limit of 10^^ erg s~^, where we only expect five or six 
background sources. The X-ray sources include one of 
the first candidate qLMXBs in an open cluster environ- 
ment (XI), a candidate CV (X2), and two active binary 
systems (X6 and X9). 

The low probability of creating a primordial qLMXB 
in a cluster the size of NGC 6819 would seem to indicate 
that source XI is most likely form ed dynamically . 

The established relationship in Pooley fc HutJ ( 2006 ) 
between X-ray sources and encounteFlrequency m glob- 
ular clusters predicts no X-ray sources in comparatively 
lower density open clusters. Even if the sources are pri- 
mordial, the known scaling of the number of primordial 
X-ray sources per unit mass in low-density globular clus- 
ters predicts no X-ray sources in the mass range of NGC 
6819. The detection of multiple X-ray sources in NGC 
6819 indicates a departure from known relationships be- 
tween X-ray sources and dynamics and mass in globular 
clusters. It may be that primordial X-ray sources are able 
to remain in an open cluster environment in comparison 
to the higher density core of globular clusters (Davies 
1997). In addition, there may be unexpected dynami- 
cal creation processes responsible for the presence of XI, 
the candidate qLMXB. We encourage further theoretical 
work in this area to help explain the creation of close X- 
ray binaries in lower density environments than has been 
previously investigated. Completion of our X-ray study 
of open clusters is necessary in order to determine the 
full extent that dynamics play in creating X-ray sources 
in open cluster environments. 
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